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Abstract 

The parametric fluorescence from a nonlinear crystal forms a conical radiation pattern. We measure the 
angular and spectral distributions of parametric fluorescence in a beta-barium borate crystal pumped by a 
405-nm diode laser employing angle-resolved imaging spectroscopy. The experimental angle-resolved spectra 
and the generation efficiency of parametric down conversion are compared with a plane-wave theoretical 
analysis. The parametric fluorescence is used as a broadband light source for the calibration of the instrument 
spectral response function in the wavelength range from 450 to 1000 nm. 
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I. INTRODUCTION 



Signal-idler photon pairs generated by spontaneous parametric down-conversion (SPDC) have 
been used to address fundamental issues of quantum theory and have found application in quantum 
entanglement and quantum information processing and metrology [IH3]. The SPDC process, also 
known as parametric fluorescence or parametric scattering [5l-U0| , is a second-order optical process 
in which a driving pump photon is scattered into signal-idler photon pairs subject to energy and 
momentum conservation. This spontaneous parametric emission can be described properly only by 
field quantization. Since its prediction and observation in the 1960s, parametric fluorescence has 
become a technique for measuring second-order nonlinear susceptibilities |10H13j and for developing 
tunable light sources via parametric oscillation or amplification processes. In 1969, Zeldovich and 
Klyshko [2] first proposed the use of parametric fluorescence (luminescence) as a nonclassical 
source of photon pairs. This description was experimentally verified by Burnham et al. in 1970 

The possible wave vectors of the signal-idler photon pairs are determined by energy and mo- 
mentum conservation, a constraint referred to as phase-matching, leading to highly directional 
parametric emission. The phase-matching condition frequently cannot be met for specific wave- 
lengths of interest or practical applications owing to the limited tunability of inherent dispersion 
of nonlinear materials. However, it can be met by selecting polarization birefringent crystals with 
appropriate refractive indices or by designing waveguides or periodic structures of specific wave- 
lengths. There are two major types of polarization phase-matching schemes for parametric down- 
conversions: Type-I, where signal-idler photons have the same polarization (co-linearly polarized 
photons), and Type-II, where the signal-idler photons have orthogonal polarization (cross-linearly 
polarized photons). Both types of parametric processes have been used to generate photon pairs, 
sometimes referred to as biphoton states, which exhibit correlation/entanglement for variables 
including polarization, momentum, time, energy, and angular momentum. 

When the phase-matching condition is met, the signal and idler radiation form a conical pat- 
tern independent of the intensity of the pump source. The angular distribution of parametric 
fluorescence is determined by the energy of the pump, signal, and idler waves, subject to the dis- 
persion of the crystal and walk-off angles of these three waves. The magnitude of the second-order 
nonlinear susceptibility is a typical selection criterion for parametric downconversion. Many 
uniaxial or biaxial nonlinear crystals have been used for parametric down conversion: for example, 
KD*P (potassium dideuterium phosphate, KD2PO4), BBO (beta-barium borate, /3 — BaB2C>4), 
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and LBO (lithium niobate, LiNbOs). In this report, we use BBO, a negative uniaxial class 3m 
crystal characterized by a wide range of transparency over the ultraviolet (A ~ 200nm) to the 
infrared (A ~ 3500nm) portion of the spectrum [16J. BBO crystals have been widely studied for 
harmonic frequency generation, optical parametric oscillation, and generation of bi-photon states. 

II. EXPERIMENTAL METHODS 

We measure the angular distribution and photon flux of parametric fluorescence from a 3-mm 
thick BBO crystal. The BBO crystal is cut at an angle of 9 m = 29 ± 0.5 ° with respect to the 
optical axis. This cut angle 6 m is chosen for the Type-I (e — > o + o) degenerate parametric down- 
conversion at A = 810 nm with a pump X p = 405 nm. The crystal is mounted on a three-axis rotary 
mount with the crystal's optical axis (OA) in the horizontal plane. The angle formed by the OA 
and the pump's propagation wave vector can be finely tuned by tilting the crystal to satisfy the 
phase-matching condition for various m near the crystal cut angle. We can thus adjust the pump 
and signal Poynting vectors from collinear to non-collinear and generate parametric fluorescence 
with varying conical emission angles. 

The pump is a violet diode laser with a TEMqq linearly polarized 2-mm 1/e 2 diameter output 
beam at a wavelength X p = 405 ram (CNI Laser MLL-III-405). The pump beam is focused on 
the crystal through a lens (LI) with a focal length of 500 mm. Lens LI and the objective are 
positioned to form a telescope such that the residual pump beam is collimated with a reduced 
beam radius below lOOum. By passing the pump beam through a pair of a half-wave plate (HWP) 
and a Glan-Taylor polarizer (PI), we can vary the incident pump intensity by rotating the HWP 
while maintaining the degree of linear polarization better than 99.9%. 

The angle-resolved images (Fig. [3]) and spectra (Fig. [6| of parametric fluorescence are measured 
by a Fourier transform optical system, including a 20 x long- working-distance objective and an 
imaging spectrometer as shown in Fig. [T] The BBO crystal is positioned at the focal plane of the 
objective lens (effective focal length f D = 10 mm). The parametric fluorescence with an amplitude 
distribution F(x, y) at the crystal is collected by a 20 x microscope objective with a 10- mm effective 
focal length (numerical aperture N.A. = 0.26). The back focal plane of the objective is the Fourier 
transform plane with coordinates (u, v) = (f a x s%Th(J)x) , fo x szra($y)). The collection angle is 
within ±15°, limited by the objective. The objective lens converges parallel rays emanating from 
the crystal to the back focal plane of the objective. In this plane, the fluorescence image in the 
crystal is transformed into a far-field image in spatial frequency that is related to the emission angle 
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as described above. The spatial intensity distribution of parametric fluorescence in the back focal 
plane of the objective lens thus corresponds to the angular distribution of radiation. This Fourier 
transform plane is placed at the front focal plane of lens L2 (focal length / = 100 mm). Lenses 
L2 and L3 are identical and separated by a distance of 2f, and they relay the Fourier transformed 
images to the entrance plane and then onto the charge-couple device (CCD) through the zero-order 
diffraction off the grating of the imaging spectrometer (PI- Acton SpectroPro 2750i, focal length 
750 mm). In this way, we measure the angular distribution of parametric fluorescence. When lens 
L2 is removed, we project the real-space spatial intensity distribution of parametric fluorescence 
from the BBO crystal onto the CCD. 

The fluorescence image is recorded by a CCD positioned in a conjugate imaging plane of the 
Fourier plane. The resultant intensity distribution is related to the Fourier transform of the inten- 
sity of the parametric fluorescence I(x,y) = \F(x,y)\ 2 . By projecting this far-field image through 
the entrance slit and the first-order diffraction of a 300 lines/mm grating, we obtain angle-resolved 
spectra as an image by taking the spectral dispersion of the parametric fluorescence as a function 
of angle. The spectral resolution of 0.1 nm is determined by the dispersion of the grating and 
the width of the entrance slit (~ 100/iin). The spatial and angular resolutions of the system are 
approximately 2 fim and 2 mrad, respectively, limited by the pixel size of the liquid-nitrogen-cooled 
CCD camera. 

In a parametric scattering process, only about one out of 10 10 incident photons is parametrically 
down-converted. It is essential to prevent the transmitted and scattered pump photons from 
entering the spectrometer. For Type-I phase matching in a negative uniaxial BBO crystal (e-o-o 
case) , the polarization of the high-frequency pump (e-wave) is orthogonal to the polarization of the 
signal and idler (o- waves). Thus, the transmitted and scattered pump photons can be suppressed 
by approximately 4-5 orders of magnitude by a pair of polarizers (PI and P2) with orthogonal 
polarization orientations for pump and single/idler waves, respectively. The pump photons are 
further rejected by a thin-film notch filter (Semrock 405-nm StopLine single-notch filter) and two 
longpass filters (a Semrock 409-nm blocking edge BrightLine long-pass filter and a Schott GG435 
glass filter). The filters are arranged in the sequence shown in Fig. [l]to suppress fluorescence from 
filters induced by the transmitted violet pump laser. The combination of polarizers and filters allows 
for a rejection of pump photons by approximately a factor of 10 10 . This rejection ratio of 10 10 can 
be further improved above 10 14 by a miniature beam blocker made of a ~ 0.5 mm-diameter silver- 
paste dot on a microscope cover positioned in front of the notch filter (NF)/objective. Depending on 
the signal wavelength, the measured angle-resolved spectra may still contain residual transmitted 
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and scattered pump and background fluorescence from filters. We measure such a background 
emission spectrum in the absence of Type-I e-o-o parametric fluorescence by rotating the BBO 
crystal 90 degrees azimuthally. 

III. EXPERIMENTAL RESULTS AND MODELING 
A. Phase Matching 

The three-wave parametric processes are calculated according to the conservation of energy 
and momentum, commonly referred to as phase matching. The angle-resolved spectra of the 
parametric fluorescence are consistent with the tuning curves calculated for the phase-matching 
condition under a plane-wave approximation. 

The energy conservation condition is expressed as 

UJ p = UJ s +UJi, (1) 

where uj p is the frequency of the incident pump wave and w s and u>i are the frequencies of the signal 
and idler waves. 

The momentum conservation condition can be expressed as 

kp = k s + kj , (2) 

where k p , k s , and kj are the pump, signal, and idler wave vectors, respectively. For Type-I down- 
conversion in a BBO, the signal and idler labels are arbitrary. In the case of degenerate down 
conversion, k s = ki, and Eq. ^ reduces to 

n p = n s cos(0' s ) , (3) 

where n p and n s are the indices of refraction of the pump and signal, and 6' s is the angle formed 
by the propagation directions of the signal and pump waves inside the crystal. 

We use a BBO crystal cut at an angle of 9 m = 29 ± 0.5°, optimized for Type-I paramet- 
ric down-conversion (e —> o + 6). Down-converted signal/idler photons are co-linearly polarized 
but orthogonal to the polarization of the pump wave. The wavelengths and wave vectors of the 
parametric fluorescence are determined by the phase-matching angle 9 m , the angle formed by the 
optical axis of the crystal (z'-axis), and the wave vector of the pump wave (z-axis) as shown in 
Fig. [2j The crystals are mounted on a rotation stage such that the optical axis lies in a horizontal 
plane when the parametric fluorescence signal is maximized. By tilting the crystal, we can vary 
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the phase-matching condition from collinear to non-collinear, leading to a conical angle up to 5 
degrees for degenerate parametric fluorescence near 810 nm. 

For Type-I phase matching, the incident pump photons are subject to the extraordinary index of 
refraction h(9 m ), while the down-conversion photons are subject to the ordinary index of refraction. 
The extraordinary index of refraction, n, depends on the phase-matching angle 9 m and follows the 
relationship: 



In the parametric process, a pump wave of wavelength X p creates signal waves at A s , and angles 
9 S , subject to energy conservation (Eq. 0) and momentum conservation (Eq. Q). In the Type-I 
e-o-o case, n s = n (X s ) and n p = n(9 m , X p ) (Eq. Q. Here the labeling of signal and idler waves is 
arbitrary (9 S = 9i). A continuum of phase-matching functions ^(A s , 9 S ) for parametric fluorescence 
can be obtained using the aforementioned equations and indices of refraction n {\ s ) and n e (X s ). 
Indices of refraction of wavelengths ranging from 0.3 /jm to 5 /im are extracted from "NIST 
Noncollinear Phase Matching in Uniaxial and Biaxial Crystals Program" as described in Ref. |17| . 
We calculate the phase matching functions (tuning curves) for down-converted signal/idler waves 
ranging from 430 to 1000 nm for a pump wave X p = 405 nm. 



B. Angle- Resolved Imaging 

We adopt the plane- wave analysis developed in Refs. |3 [181 EH] to determine the angular 
distribution of parametric fluorescence. The effects of a finite pump beam size have also been 
considered in, for example, Refs. \20\ 12 lj . Under a plane- wave approximation, the parametric 
fluorescence forms a conical angular distribution The diameter and axis of the conical emission are 
determined by the wavelengths of the pump, signal, and idler waves, the Poynting vector walk- 
off angles of the three interacting waves, and the dispersion of the nonlinear crystal. In such a 
parametric process, the angular spread of the conical emission is determined by the conservation 
of transverse and longitudinal momenta of interaction waves. The transverse momentum induced 
by focusing the pump into the crystal contributes to a finite angular spread of the cones. By 
considering the longitudinal momentum, we deduce that parametric fluorescence signal is inversely 
proportional to the interaction length. In our experiments, we use a lens with a long focal length 
(/ = 500 mm), resulting in a focal spot with a 1/e 2 radius larger than 100/im. Thus, for the 
experiments reported here, the effects of walk-off and finite pump beam size are negligible compared 
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with the spectral and angular resolution of the optical system. 

We can determine the angular distribution of parametric fluorescence using the following phase- 
matching function for a finite crystal length L and a pump Gaussian beam profile with a 1/e 2 radius 
W H3E21 [23]: 




The mismatch wave vector, Ak, is decomposed into longitudinal (z || k p ) and transverse (in 
x-y plane) parts: Ak z and k = \J Ak x + Ak y . The phase-matching tolerances can be considered 
in terms of the angular spread (A9 S ) and spectral bandwidth (AA S ), defined as the full-width-at- 
half-maximum (FWHM) for the above function. For a pump wave with a focal radius W ~ 50/im, 
only the tolerances from the sine 2 part can be measured in our optical system. For this situation, 
considering a Taylor series expansion of Ak z near the perfect phase-matching point (Ak z = 0), we 
can analytically deduce the angular spread (A#fwhm) and spectral bandwidth (AAfwhm) for the 
degenerate case (k s = ki): 



2 x 0.886^ 2.783 x n s 

A0 FWH m(0 s ) = L ~ ]dAkz/d9sl ~ L x k s x 9 S ' (6) 

and 

2 x 0.886vr 0.443 x A s 2 

s Lx\dAk z /d\ s \ L x n s cos(8 s /n s ) ' 

where Ak z = \k p — 2k s cos(# s )|, kj = n s uj s /c = 2-7rn s /A s , and 9 S <^ 1. 

The parametric fluorescence signal of angle-resolved images at A = 810 nm are shown in Fig. [3} 
These false-color images, taken through a 1-nm band-pass filter, represent the angular intensity 
distributions of parametric fluorescence at A = 810 ± 0.5 nm for the phase-matching angle m — 
28.6°, 28.8°, 29.1°, and 29.4°. The BBO crystal is cut at the designed angle with about 1° tolerance. 
To determine the phase-matching angle 9 m with better precision, we first set the phase-matching 
angle for the collinear case by comparing the simulated and experimental angular distributions. 
We then deduce the phase-matching angle for the non-collinear case from the tilting angle of the 
crystal relative to that for the collinear case. The conical signal angle (9 S ) of degenerate parametric 
fluorescence at A = 810 nm increases with 9 m . In Fig. El we plot the angular spread (A#fwmh) 
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as a function of the inverse of the signal angle (1/9 S ). The angular spread decreases with 9 S for 
small angles when sin(# s ) ~ 9 S . We attribute the discrepancy between experimental data and Eq. 
([6]) to a limited experimental angular resolution (~ 2 mrad), finite pump beam size and spatial 
coherence, and the birefrigent walk-off. 

C. Angle-Resolved Spectroscopy 

The parametric fluorescence flux per unit frequency is 

N s (u s , k.) = n d l f ( f Ui Up 12 Np du> s d 2 n s [ exp( V) • sine 2 (\ L Ak z ) , (7) 
2Tr*c 6 e n s ninp J 2 \2 J 

where d e jf is the effective second-order nonlinear coefficient, L the interacting crystal length, 
N p = Pp/hujp the pump flux, £ = kW = (n s + Kj) W the dimensionless transverse momentum 
associated with the signal and idler waves, and W the 1/e 2 pump beam radius. The integration 
over £ is detailed in Appendix. 

Experimental angle resolved spectra N*(X S ,9 S ) and corresponding calculated N S (X S ,9 S ) are 
shown in Fig. [6](e)-(h). The tuning curves, corresponding to the perfect phase- matching condition, 
are shown as white dashed lines on the experimental spectra. The experimental parameters are 
N p = 1.63 x 10 17 /s (P p = 80 mW ), L=3 mm, and W 2 = 60^m x 30^m. The indices of 
refraction n s , and n p are evaluated using the database in Ref. |17j (n Sj j ~ 1.66 at 810 nm). 
The effective nonlinear coefficient of a BBO crystal can be deduced from its d-matrix using d e ff = 
d^i s'm(9 m + p) — d22 cos(# m + p) sin(3^>) for Type-I phase matching. 9 m is the phase-matching 
angle, <j) the azimuthal angle, and p the birefringent walk-off angle. We use d e ff =1.75 pm/V for 
405 nm adopting from Refs. [U IT3tl2"l"]. 

Computationally, the angular distributions of parametric fluorescence are calculated for signal 
wavelengths from 430 to 1000 nm with step size 5X = 1 nm and 59 s = 0.03° ~ 0.5 mrad. Selected 
calculated angular distributions of fluorescence flux per 1 nm are shown in Fig. [6] (e)-(h) for 9 m = 
28.6°, 28.8°, 29.1°, 29.4°. The shortest signal wavelength appears in the simulation is ~ 440 nm, 
limited by the availability of refractive indices between A = 0.3 pm and 5 pm [25J. Experimentally, 
parametric fluorescence with a wavelength as short as 431 nm can be observed near 9 S = 0. 

The experimental angle-resolved parametric fluorescence images are shown in the left panel of 
Fig. [HJ These images are taken through a notch filter and longpass filters (cutoff wavelength of 420 
nm) as shown in Fig. IT} Experimental angle-resolved spectra are acquired by spectrally resolving 
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the parametric fluorescence across the fluorescence cone center through the entrance slit with an 
opening of 100 /im, corresponding to A9 X ~ 0.6°. Angle-resolved fluorescence images (Fig. [6](a)- 
(d), left panel) are taken for the zero-order diffraction of a holographic grating with 1800 lines/mm, 
while angle-resolved spectra (right panel) are dispersed by a grating with 300 lines/mm and a blaze 
wavelength of 1000 nm. The angular resolution is about 2 mrad, while the spectral resolution is 
about 0.1 nm. Residual stray or scattered pump laser light can be measured by rotating the BBO 
crystal 90° azimuthally. Such background 'noise' is subtracted. Note that the measured CCD 
intensity is subject to the nonuniform spectral sensitivity and collection efficiency of the optical 
system including the effects of lens coating, optical filters, and gratings, and the spectral response of 
the CCD camera. The spectra branches out from 6 S = at A ~ 433 nm. The secondary weak arcs 
in the wavelength range above 870 nm outside of the expected parametric fluorescence peaks are 
due to the second-order diffraction of the grating for the fluorescence from approximately 435 nm 
to 500 nm. The inner arcs branching from 435 nm and closing near 530 nm are due to parametric 
fluorescence from a Type-II parametric process (e — > o + e). The turning curves for such Type-II 
parametric fluorescence are indicated by black dashed-doted lines in experimental spectral images. 



D. Fluorescence Photon Flux 



The integrated parametric fluorescence photon flux can be obtained by the integrating over £ 
and k s in Eq. ^ . For values of 6 m or co s such that the parametric fluorescence cone has a radius 
sufficiently large with negligible emission at the cone center, the resultant integrated fluorescence 
flux is p2] 

it eo n p z 
= (27r)4 McdgL 

The efficiency rj s = N s /N p is a coefficient depending largely on the material properties such as the 
second-order nonlinear coefficient, interacting crystal length, and index of refraction for the pump 
wave. Assuming a bandwidth AA = 1 nm and L = 3 mm, the efficiency coefficient rj s = 1.3 x 10~ 10 
for the degenerate parametric fluorescence at X s = 810 nm under X p = 405 nm. Specifically, we 
evaluate r] s for 9 m = 29.12° (Fig. [6^). We integrate the photon flux of simulated angle-resolved 
spectra for X s = 809.5 — > 810.5 nm, 9 S = — > 7.5°, and eft = — > 2tt. The total parametric 
fluorescence photon flux, including both degenerate signal and idler waves, is 2N S ~ 4.2 x 10 7 /s . 



The wavelength of parametric fluorescence generated here ranges from ~ 430 11m to above 1000 
nm. The collection efficiency and spectral response of the optical systems could vary more than 
an order of magnitude in such a broad wavelength range. Therefore we consider the degenerate 
parametric fluorescence at A = 810 nm to compare the experimentally determined fluorescence 
photon flux with the theoretically integrated photon flux (see Eq. Q). 

The integrated degenerate parametric fluorescence flux at A s = 810 nm as a function of the 
incident pump flux are shown in Fig. IU First, we determine the overall collection efficiency and 
the spectral response of the imaging and spectroscopy system by passing a laser beam of A = 810 
nm with a known photon flux through the same optical path as the parametric fluorescence. We 
then integrate the signal in angle-resolved images as shown in Fig. [3j Taking into account the 
system response at A = 810 nm, we can determine the absolute fluorescence flux experimentally 
within roughly 20% error. The relative fluorescence fluxes, determined by the linearity of the CCD 
camera, is within 1% error. The absolute value of pump flux is within roughly 10% error. The 
relative pump fluxes, as varied by a combination of a half- wave plate and a polarizer and limited by 
the linearity of the power meter, are known within a few percent. Thus we can investigate the pump 
flux (power) dependence of parametric fluorescence flux with precision. Fluorescence signal flux is 
linearly proportional to the pump flux over two order of magnitude, confirming that the dominant 
signal is spontaneous parametric fluorescence as described by Eq. Q. The slopes correspond to 
the efficiency coefficients. Considering both signal and idler fluorescence, we determined r\ = 2r) s 
to be approximately 2.8 x 10~ 10 . Experimental and theoretical values of n = 2n s for selected 
phase-matching angles are listed in Table |I} 

E. Instrument Spectral Response Function 

Parametric fluorescence spectra can also be used to calibrate the instrument spectral response 
function (ISRF) of the imaging spectroscopy system. According to Eq. ([8]), which is valid for non- 
collinear cases, we can deduce a parameter S = N s x Af A^ = 2 (27r) 4 U cd e ff 2 L N p dX s /(eo n p 2 ) 
|18j . S is a wavelength-independent constant for a given pump wavelength and geometry. We 
define a generalized spectral function, S(X S ) = N S (X S ) A 4 A?, for both calculated and experimental 
angle- resolved spectra. The calculated S s i m (X s ) exhibit less than 1% variation between A = 460 
and 1000 nm. Experimentally, S* (A s ) = N*(X) A^Af can be determined from the integration of 
an angle-resolved spectrum N*(X, 9 S ) over 6 S . S exp (X) represents the relative ISRF of the imaging 
spectroscopy system, including the optical components, grating, and CCD camera along the fluo- 
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rescence collection optical path. The value of the ISRF at a fixed wavelength can then be used to 
determine the absolute ISRF across the parametric fluorescence wavelength range. The effective 
efficiency, defined as (# of photo-generated electrons / # of photons), is «20% at A = 810 nm in 
our experiments. S* xp (X) and S S i m (\) for phase-matching angles 9 m = 29.1° and 29.4° are shown 
in Fig. [7J The stray scattered pump laser signal becomes increasingly difficult to subtract from 
these angle-resolved spectra, leading to a distorted S*(X). We thus use S* xp (\) at 6 m = 29.1° to 
deduce the ISRF of our optical setup shown in Fig. [TJ 
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IV. APPENDIX: CALCULATION OF FLUORESCENCE FLUX 

Here we describe the integration over £ of Eq. ([7]) for the calculation of the angular distribution 
of fluorescence flux here. The mismatch wave vector Ak = k' s + k\ — k p can be decomposed into 
longitudinal (Ak z z) and transverse {k) parts: 

Ak z = k sz + k iz — k pz 



K 2 -Ks + \/kf-rf-k' pz , and 

Here k' s = n s uj s /c, k\ = riiUi/c, and k' p = n p uj p /c are the wave numbers for the signal, idler, 
and pump, respectively. k s (k^) is the transverse wave vector of the signal (idler) wave. The 
phase- matching function $ in Eq. ^ is a function of three variables: uj s , k s and Kj. Considering 
energy conservation uj p = u s + cjj, we can carry out the integration over £ for two independent 
variables, oj s and k s . Using £ = £ s + £j = (k s + Ki) W, we rewrite Eq. Q as 



AT I \ Hde ff U s UJiUJpL 2 N p 2 f 2 / 1 2 \ . 2 /l \ 

N s {u) a ,K s ) = JJ v - (hj s d k s I d ^ exp -- £ s + & -sine -LAk z 

87T 4 c s e n s in n p J \ 2 J \2 J 

We further simplify the numerical integration by (a) considering that the sine 2 term is a constant 
for a given set of uj s and k s , and (b) applying a saddle-point approximation. 
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k s and K{ form an angle <p in the x — y plane of the crystal, where <p = in the anti-parallel case. 
Under a pump wave vector along z, the phase-matching condition is met mostly for k' s + k\ = 0; 
i.e., ip ~ 0. The Gaussian term of the integrand can thus be separated and integrated with a 
saddle-point approximation for ip rj and £ s « 



i ^ 2 



y d 2 & exp (-^(£ 2 + e 2 -266 cos(p))) « J exp f-i (6 - £i) 2 J / ^exp £ 

The integrand is a maximum at £ s =6) where pairs of signal and idler photons are emitted in 
approximately opposite conical directions. Applying the above two approximations, we obtain 



N.(u a , k s , 4>) = e J\ 3 p p -V2Ttduj s k s cLk s d4> / d£i e-3^-w • sine 2 - L Ak z , 

87r 4 c 6 e n s ni n p 7 V 2 / 

where d dn s d<p, 4> is the azimuthal angle. 

N s (cj s , k s , 0) is isotropic in 0. It can be expressed as N S (X S ,9 S ), a function of experimentally 

measurable signal angle 9 S and wavelength A s by considering that k s = k s sin^ and oj s = 2ttc/X s : 



N s (X s ,e s ) = ^^Kd ef /u L*N P sin{2es)dXsd0g I . sinc2 ( i LAk A . (9) 

The equation for N S (X S ,6 S ) above, together with u p = u s + uji and k p = k s + kj, is used in 
the numerical integration to obtain the theoretical angular distribution of parametric fluorescence 
shown in Fig. [6] 
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TABLE I: Parametric fluorescence efficiency 77 = 2N S /N P for 810 ± 0.5 nm 





Exp. 


Sim. [Eq. (9|] 


Th. [Eq. (8)] 




28.6° 


1.3 x 10~ 10 


2.5 x 10~ u 




1.6612 


28.8° 


2.8 x 10~ 10 


2.5 x 10~ 10 


2.6 x lO" 10 


1.6608 


29.1° 


2.7 x 10~ 10 


2.6 x 10~ 10 


1.6602 


29.4° 


2.7 x 10~ 10 


1.6599 
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FIG. 1: Schematic diagram of the experimental setup. The transmitted and scattered pump photons 
are rejected by a miniature beam blocker, a thin-film notch filter, and long-pass filters. The angular and 
spectral distributions of conical parametric fluorescence are measured by an imaging spectrometer through 
a Fourier transform optical system. LI, L3 and L4 are convergent lenses with focal length f = 500, 100, 
and 100 mm, respectively; L2 is an objective (20x N. A. =0.26) with an effective focal length of 10 mm 
(Mitutoyo Plan Apo infinity-corrected long-working-distance objective); PI and P2 are Glan-Taylor and 
Glan-Thompson polarizers; Ml and M2 are silver mirrors; HWP is a half-wave plate for A = 405 nm; BB 
is a miniature pump beam blocker; NF is a notch filter (Semorck 405-nm StopLine single-notch filter); and 
LP represents longpass filters (a Semrock 409-nm blocking edge BrightLinc long-pass filter and a Schott 
GG435 glass filter). When L2 is removed, a real-space fluorescence image is formed at the entrance of the 
spectrometer with an imaging magnification of 10 x. Examples of real-space and angle- resolved fluorescence 
images are shown with actual dimensions. 
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FIG. 2: Schematic diagram of the crystal and laboratory frame coordinates for parametric 
down-conversion in a BBO crystal. The phase-matching angle 8 m is defined as the angle formed by the 
crystal optical axis (V) and the pump wave vector (z). The angles 9' s and 9 S are, respectively, internal and 
external angles formed by the signal and pump wave vectors. Here the incident pump wave is horizontally 
polarized, leading to a vertically polarized down-converted signal and idler waves. 
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FIG. 3: Angular intensity distribution of parametric fluorescence at 810 nm. (a)-(d) Angle- 
resolved images of parametric fluorescence for A = 810 ± 0.5 nm and 9 m as indicated. The external angle 
9 S is formed by the pump and signal wave vectors in air (see also Fig. [2]). The color palette represents 
the intensity of the single radiation. Parametric fluorescence is spectrally filtered through a 1-nm bandpass 
filter with central wavelength A = 800 nm. The phase-matching angle is adjusted by tilting the BBO crystal 
with respect to the pump wave vector. The collinear phase-matching angle is set to 6 m — 28.6° according 
to a theoretical calculation using indices of refraction given in Ref. [17]. (e)-(g) Experimental (blue dashed 
line) and theoretical (red solid line) cross-sections. The non-collinear phase-matching angle relative to the 
collinear one can be determined experimentally by measuring the tilting angle of the crystal surface with 
respect to the pump wave vector. 
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1/0 S (rad 1 ) 

FIG. 4: Angular spreads A6*fwhm of parametric fluorescence at 810 nm. A6>fwhm is plotted as a 
function of the inverse of the signal angle (1/6 S ). Experimental data are represented with error bars as solid 
circles. The red solid line is the theoretical curve according to Eq. |6j while the dashed line is the theoretical 
curve including a finite angular resolution of 2 mrad. Selected experimental angular intensity profiles for A 
(6 S = 1.3° = 0.023 rad), B (2.8° = 0.049 rad), and C (3.7° = 0.065 rad) are shown in the inset. 
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FIG. 5: Parametric fluorescence flux at 810nm. Integrated degenerate parametric fluorescence flux 
at X s — 810 nm as a function of the incident pump flux for the collinear case and three phase-matching 
angles 9 m of the angle-resolved images shown in Fig. [3j Signal flux is linearly proportional to the pump flux 
over two order of magnitude, confirming that the dominant signal is spontaneous parametric fluorescence as 
described by Eq. Q. The slopes rj = N s /N p are 1.16 x 1CT 10 for 9 m = 28.6°, 2.6 x 1CT 10 for 6 m = 28.8°, 
and 2.8 x 1(T 10 for 9 m = 29.1° and 29.4°. 
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FIG. 6: Angle-resolved spectra of parametric fluorescence, (a)-(d) Experimental angle-resolved 
parametric fluorescence images (left panel) and spectra (right panel) for a phase-matching angle 9 m = 
28.6°, 28.8°, 29.1°, and 29.4°, respectively, (e)-(h) Theoretical fluorescence flux calculated according to 
Eq. ([7]) (see also Appendix). The tuning curves for the perfect phase-matching condition are indicated by 
the white dashed lines on experimental imaging spectra. The color palette represents the calculated photon 
flux with AA S = I nm and A9 S = 0.5 mrad on a logarithmic scale. 
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FIG. 7: Instrument spectral response function. The normalized spectral density function S(X) for 
6 m = 29.1° and 29.4°. S S i m is determined from the integration over 9 S = —7.5 to 7.5° of the calculated angle- 
resolved spectra as shown in Fig.[6](Eq. The black and red curves are S s i m or S* xp for 6 m — 29.4° and 
29.1°, respectively. The value of S S i m is a constant with 1% standard deviation across wavelengths from 500 
to 1000 nm, validating the calculated angular fluorescence spectra and Eq. (JsJ) . S* xp (X) = N*(X) A s 4 Aj 2 , 
where N*(X) is the integration over 8 S « —15 to 15° of the experimental imaging spectra N* (X S ,9 S ). 
S* xp (X) represents a relative instrument spectral response function (ISRF) of the optical spectroscopy system, 
including optical filters and a Glan-Thompson polarizer along the path of the fluorescence , a liquid- nitrogen- 
cooled CCD (PI-Acton Spec-10:400BR), and a 300g/mm plane ruled reflectance grating with 1000-nm 
blaze wavelength (PI-Acton 750-1-030-1). The polarization of the fluorescence is vertically polarized. The 
absolution ISRF (or collection efficiency) can be deduced by calibrating S* xp (X) at a fixed wavelength such 
as A = 810 nm in our experiments. 
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